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We report first measurements of e + e~ —* £)M+£)(*) - processes far above threshold. The cross- 
sections for e + e~ — > Dp~ D*^~ and e + e~ — + D + D^T at %/s = 10.58 GeV/c 2 are measured to be 
0.55 ± 0.03 ± 0.05 pb and 0.62 ± 0.03 ± 0.06 pb, respectively. We set upper limits on the cross- 
sections for e + e~ — > D*J~ T)*7 , e + e~ — > D*^ D*^~ , e + e~ — > D + D*^~ and e + e~ — -> D + D~ processes. 
The analysis is based on 88.9 fb _1 of data collected by the Belle experiment at the KEKB e + e~ 
asymmetric collider. 
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The processes e + e~ — > D^*'D^*\ with no extra fragmen- 
tation particles in the final state, have not previously 
been measuered for energies ^/s » 2Mb- The cross- 
sections for these processes can be computed once the 
charmed meson form factors are determined for the ap- 
propriate value of momentum transfer, q 2 = s. In the 
HQET approach based on the heavy-quark spin sym- 
metry, the heavy meson form factors can be expressed 
in terms of a universal form factor, called the Isgur- 
Wise function. However, for large q 2 , the leading- twist 
contribution, which violates heavy-quark spin symme- 
try, becomes dominant Q. For an intermediate range of 
momentum transfer, the Isgur-Wise contribution (hav- 
ing subleading twist) is also important. A calculation 
that takes these effects into account Q predicts that the 
cross sections for e+e~ — > DD* and e + e~ — > D* L D"^ 
are equal to each other, and are about 2.5 pb at y/s ~ 
10.6 GeV (the subscripts indicate longitudinal [L] and 
transverse [T] polarizations of the £)*). The cross-section 
for e + e~ — > DD is expected to be suppressed by a factor 
of ~ 10~ 3 . The cross-section predictions are subject to 
large uncertainties (~ a factor of 3), because the theoreti- 
cal formulae contain many poorly known parameters; the 
prediction that the DD* and D* L D™ cross-sections should 
be equal is expected to be robust Q ■ Recently, new calcu- 
lations in the framework of the constituent quark model 
have become available ,3j . 

In this paper, we present the first measurement of the 
e+e~ — * D* + D*~ and e + e~ — * D + D*~ cross-sections 
and polarizations at s/s ~ 10.6 GeV. We also set an 
upper limit on the cross-section for e + e~ — * D + D~. 
The present study is limited to final states that contain 
charged D^ mesons only. Since the contribution of the 
electromagnetic current coupled to light quarks is negli- 



gible compared to that for heavy quarks, the neutral and 
charged charm meson cross-sections are expected to be 
the same Q. 

The analysis is based on 88.9 fb _1 of data at the 
T(4S f ) resonance and nearby continuum, collected with 
the Belle detector Q at the KEKB asymmetric energy 
e + e~collider We select well-reconstructed tracks con- 
sistent with originating from the interaction region as 
charged pion candidates. Those passing particle iden- 
tification cuts based on dE/dx, aerogel Cerenkov, and 
time-of-flight information Q are selected as charged kaon 
candidates. We then reconstruct D° and D + mesons 
in the decay modes D° — > K~tt + , D° — ► K~ir + ir + ir~ 
and D + — > K~tc + tt + . The selected combinations are 
constrained to a common vertex, and quality cuts are 
imposed on the vertex fit to reduce the combinatorial 
background. A 15MeV/c 2 interval around the nominal 
D masses is used to select D° — > K~ir + and D + — > 
K~7T + ir + candidates; for the D° — > K~ir + 7r + ir~ de- 
cay mode the signal window is chosen to be 10MeV/c 2 
around the nominal D° mass (~ 2 a each case) . The se- 
lected D candidates are then subjected to a mass and 
vertex constrained fit to improve their momentum reso- 
lution. The D* + mesons are reconstructed in the D°ir + 
decay mode. The mass of the D°ir + combination is re- 
quired to be within a 2MeV/c 2 (~ 3er) mass interval 
around the nominal D* + mass. 

The processes e + e~ — > £)W+DW" can be identi- 
fied by energy-momentum balance in fully reconstructed 
events that contain only a pair of charm mesons. How- 
ever, the small charm meson reconstruction efficiency of 
the studied channels results in a tiny total efficiency in 
this case. Because of the simple two-body kinematics, it 
is sufficient to reconstruct only one of the two charmed 
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mesons in the event to identify the processes of inter- 
est. For simplicity, we refer to the fully reconstructed 
D meson as the £)W + , and the other as the D*~; the 
charge-conjugate modes are included in the analysis. We 
choose the mass of the system recoiling against the re- 
constructed Of rccoi i(L> w +)) as a discriminating 

variable: M rccoi i(£> ( * )+ ) = ^ {s/s - E DM + ) 2 - pg M+ , 

where y/s is the total center of mass (CM) energy, and 
E D (,)+ and p D (,)+ are the CM energy and momentum 
of the reconstructed £)(*)+ . For the signal, a peak 
in the M rcco ii distribution around the nominal D~ or 
D*~ mass is expected. This method provides a sig- 
nificantly increased efficiency, but also a higher back- 
ground, in comparison with full event reconstruction. For 
the e+e~ — > D + D*~ and e+e~ — > D* + D*~ processes 
we find a better compromise between higher statistics 
and smaller background: the first D^*> + is fully recon- 
structed, while the recoiling D*~ is required to decay into 
D°ir~ low . The reconstructed tt~ 1ow provides extra infor- 
mation that allows us to reduce the background to a neg- 
ligible level using the difference between the masses of the 
systems recoiling against the D^ + tt~ 1ow combination, 
and against the £)(*)+ alone, AM rocoil = M recoil (£) ( * )+ )- 
M Teco i\(D(*} + Tr~ low ). The variable AM recoil peaks around 
the nominal D*~—D° mass difference with a resolution 
of cr A M„ coil ~ 1 MeV/c 2 as found by Monte Carlo (MC) 
simulation. For e+e~ — > D* + D*~ and e+e~ — > D + D*~ 
we require AM reC oii to be within a ±2MeV/c 2 interval 
around the nominal M D *- — M-go mass difference. 

The M Teco [i(D* + ) and M reco ii(-D + ) distributions are 
shown in Figs.^, and^ 3 , respectively. Clear signals are 
seen around the nominal D*~ mass in both cases. The 
higher recoil mass tails in the signal distributions are due 
to initial state radiation (ISR). The M recoil distributions 
for events in the AM rocoil sideband (0.150 MeV/c 2 < 
AM rcco ii < 0.154 MeV/c 2 ) are shown as the hatched his- 
togram (barely visible in Fig.Q^ due to its small size). 

The backgrounds in the region M reco n < 2.1GeV/c 2 
are negligible for both processes, so we consider this in- 
terval as the signal region. There are three possible back- 
ground sources: 

I incorrectly reconstructed D* + or D + ; 

II e + e~ — > £>(*)+ Dnir~ (n > 0), where the tt~ 1ow is 
not produced from D*~ decay (and can be either 
from fragmentation or from unreconstructed D de- 
cay), and thus produces no peak in the AM reco ii 
distribution; 

III e+er -> DW+D*-mr, where n > 1. 

First we consider the process e + e~ — > D* + D*~. To 
estimate background (I) numerically, we count the en- 
tries in the signal region for D°ir + combinations taken 
from the D*+ mass sideband (2.016 GeV/c 2 < M D o w+ < 
2.020 GeV/c 2 ). Three events are found in the data, while 
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FIG. 1: a) M Icco n(D* + ) and b) Mr eC oii(-D + ) after applying 
the AM reco ii requirement. Points show the AM rC coii signal 
region; hatched histograms show the AM re coii sideband; solid 
lines represent the fits described in the text; dashed lines show 
the contribution due to events with ISR photons of significant 
energy; dotted lines are the expected background contribu- 
tion, c) M reco ii (D + ) after D + sideband subtraction without 
requiring an extra ttJ 1ovj in the event. 



the MC predicts a contribution of 2.5 events from the 
signal process due to non-Gaussian tails in the M D ,+ 
resolution function. The signal MC is normalized to 
the number of entries in the M rcco ii(-D* + ) < 2.1 GeV/c 2 
region in the data. Background (II) is estimated us- 
ing the AM reC oii sideband (0.150 MeV/c 2 < AM rocoil < 
0.154 MeV/c 2 ). In the signal region 8 events are found; 

4 events are expected according to MC from the signal 
process because of significant energy ISR. Thus back- 
grounds (I) and (II) are estimated to be smaller than 

5 and 10 events at the 90% CL, respectively. The re- 
maining background (III) can result in peaks in both the 
M{D* + ) and AM rcco ii distributions, but has a threshold 
in the M reco ii distribution at M D *++M„a = 2.15 GeV/c 2 , 
which is ~ 1 a higher than the chosen M reC oii signal in- 
terval. To estimate the residual background (III) con- 
tribution in the signal region we perform a fit to the 
^recoil (D* + ) distribution. The signal function is de- 
termined from the MC simulation and parametrized as 
the sum of a core Gaussian and an asymmetric func- 
tion representing the case when the studied process is 
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accompanied by radiative photon(s) with significant en- 
ergy (E IS r > lOMeV). The Al rccoi i(L>* + ) resolution due 
to detector smearing and the signal function offset are left 
as free parameters in the fit to check the agreement with 
the MC predictions. The background (III) distribution is 
parameterized by a threshold function, a(M roco ii(Z?* + ) — 
M(D*-) PDG - M(ir°)pDG) , convolved with the detec- 
tor resolution, where a and (3 are free parameters. The 
fit results are shown in the Fig. ^ as a solid curve; the 
dashed line shows the contribution of the studied process 
with significant energy ISR photons (Ejsr > lOMeV). 
The dotted line represents the expected background (III) 
distribution. The signal yield is found to be 815 ± 28 
events in the whole fit region (676 ± 24 events in the sig- 
nal region M reC oii < 2.1 GeV/c 2 ). The Af rcC oii resolution 
a = 56.1 ± 2.2MeV/c 2 is found to be in excellent agree- 
ment with the MC expectation (56.4 MeV/c 2 ), and the 
shift of the signal peak position in the data with respect 
to the MC position is found to be consistent with zero 
(0.6 ± 2.5MeV/c 2 ). The contribution from background 
(III) in the signal region is estimated from this fit to be 
less than 2 events at the 90% CL. 

For the process e + e~ — > D + D*~ , proceeding in a sim- 
ilar way, we find a signal yield of 423 ± 20 events in the 
whole fit region (360 ± 17 events in the signal region 
Mccoii < 2.1GeV/c 2 ), with backgrounds (I), (II) and 
(III) smaller than 8, 6 and 2 events at the 90% CL, respec- 
tively. Finally we estimate the total background in the 
M reco n < 2.1 GeV/c 2 interval to be smaller than 13 and 
10 events for the e+e~ ->■ D*+D*- and e+e" -> D+D*- 
processes, respectively, which is of the order of 1% of the 
signal. We therefore assume that all events in the interval 
A/rocoii < 2.1 GeV/c 2 are signal, and include the possible 
background contribution in the systematic error. 

Since the reconstruction efficiency depends on the pro- 
duction and D*^ helicity angle distributions, we per- 
form an angular analysis before computing cross-sections. 
The helicity angle of the non-reconstructed D*~ is cal- 
culated assuming two-body kinematics. A scatter plot of 
the helicity angles for the two D*-mesons from e + e _ — > 
D*+D*- (cos(f>{D* rec ) vs. cos<f>(D* on _ rec )) for the sig- 
nal region is shown in Fig. E^. This two dimensional 
distribution is fitted by a sum of three functions corre- 
sponding to the D^D^, D* T D* L and D* L D* L final states, 
obtained from the MC simulation. The fit finds 6^3, 
708 ± 36 and 4±^f events associated with D* T B* T , D* T D* L 
and D* L D* L final states, respectively. Figured shows the 
D*~ meson helicity distribution for e + e~ — > D + D*~. 
A fit finds 433 ± 24 and — 1.5 ± 2 events corresponding 
to DDIj, and DD^, respectively. We conclude that in 
e+e~ -> D*+D*- and e+e~ -> D+D*~ production the 
final states are saturated by D^D* L and DD^. (DD^ is 
required by angular momentum and parity conservation 
for e + e~ — * D + D*~ via a virtual photon; the D^D* L re- 
sult is non-trivial.) The production angle distributions 
for e + e~ — > DX^D* r ~ and e + e~ — > D + DXr are, there- 
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FIG. 2: a) A scatter plot of cos(0o» ec ) vs. cos((f>o* _ ) 
for e + e~ — > D* + D*~events. b) D*~~ meson helicity angle 
distribution for e + e~ — > D + D*~ signal candidates. The curve 
represents the fit described in the text. 



fore, fixed to be 1 + cos 2 9 in both cases. As a cross-check 
we study the production angle distributions for D* + from 
e+er -> D*+D*- and D+ from e+e" -> D+D*~ pro- 
cesses. After correction for reconstruction efficiency in 
bins of the production angle we fit the production angle 
distributions with the function (1 + Acos 2 6>). The pa- 
rameters A are found to be equal to 0.79^0 30 anc ^ 2.3^q'y 
for the two processes, which are in agreement with the 
expected value A = 1 for both processes. 

To calculate the Born cross-section for the studied pro- 
cesses we determine the fraction of events in the signal 
region with an ISR photon energy smaller than the cho- 
sen cutoff (.©cutoff = lOMeV) using a MC simulation. 
In the MC we assume a 1/q 6 dependence of the ratio of 
the cross section to er(e + e~ — > /« + M~) as predicted by 
Ref. 0. We also try a 1/g 2 dependence (corresponding 
to flat form-factors) and include the resulting shift of the 
event fraction (< 1.3%) in the systematic error. The 
reconstruction efficiencies are determined from MC sim- 
ulation. The Born cross-sections for e + e~ — > £)(*)+D*~ 
are calculated according to the formulae in Ref. and 
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are listed in Tabled The final result is independent of 
the choice of E cuto e- The sources of systematic error 
are summarized in Table [n] The dominant contributions 
are from the uncertainties in tracking efficiency and D^*> 
branching ratios. 



Process 


OBorn (pb) 


e 


~e~ 


-> D* T 




< 0.02 @ 90% CL 




~e~ 


-> D* T 


¥d T 


0.55 ±0.03 ±0.05 


e- 


~e~ 


^Dt 


'D* L - 


< 0.02 @ 90% CL 


<r 


~e~ 


^D+ 


Df 


< 0.006® 90% CL 




~e~ 


-^D+ 




0.62 ±0.03 ±0.06 




~e~ 


-^D+ 


D~ 


< 0.04 © 90% CL 



TABLE I: The Born cross-section results. The first error is 
statistical and the second one is systematical. 



TABLE II: Summary of systematic errors for the e+e 
D* + D*~ and e+e~ — > D + D*~ measurements. 



Source 


e + e" -> D* + D*- 


e + e" D + D*- 


Tracking efficiency 


7% 


5% 


Identification 


2% 


2% 


Backgrounds 


+0 % 

-1.6 70 


t\. 5 % 


Form-factor shape 


1.3% 


1.3% 


Luminosity 


1% 


1% 


B{D M ) 


4% 


8% 


Total 


9% 


10% 



We search for the process e + e~ — * D+D~ by studying 
the mass spectrum of the system recoiling against the re- 
constructed D+ without requiring an extra soft pion in 
the event. In the e+e" -> D*+D*~ and e+e" D+D*~ 
analyses, backgrounds are strongly suppressed by the 
AM roco ii cut, which is not applicable for the e+e" — > 
D+D~ search; without this requirement the combinato- 
rial non-D + background is significant. We use D+ mass 
sidebands (20MeV/c 2 < \M K ^ - M D \ < 35MeV/c 2 ) 
to extract the M rcco ii distribution for the combinato- 
rial background. Figure lc shows the M Teco i\(D+) dis- 
tribution after D+ mass sideband subtraction. We 
fit this distribution with the sum of two signal func- 
tions corresponding to D~ and D*~ peaks and a back- 
ground function. The latter is a threshold function, 
a(x — M(D~)pdg — M( 7T °)pdg) 13 , convolved with the 
detector resolution, where a and (3 are free parameters. 
For the fit we use only the region M rccoi i < 2.25 GeV/c 2 , 
because of a possible contribution of e+er — > £)(*)£)** at 
higher M re coil- The fit finds 13±24 events in the D~ peak 
and 935±42 in the D*~ peak. The fit function is shown in 
Fig. lc as the solid line; the dashed line shows the contri- 
bution of events with ISR photons of significant energy 
(larger in this case due to the absence of the AM reco ii 
cut); and the dotted line represents the case where the 
contribution of e+e~ — ► D+D~ is set at the value corre- 
sponding to the 90% CL upper limit. The reconstruction 



efficiencies for e + e~ — * D+D~ and e+e~ — ► D + D*~ 
are found from MC. The production angle distribution 
for e + e~ — * D+D~ is assumed to be proportional to 
sin 2 9, while the production angle for e+e~ — ► D+D*~ 
is fixed from the study with the AM rGCO ;i requirement. 
The e + e~ — > D+D*~ Born cross-section is calculated to 
be 0.54 ± 0.04 pb. In this method the systematic uncer- 
tainty in the signal yield is larger than in the AM reco ii 
method due to significant e + e~ — > D+Dir background 
under the peak, which can only be extrapolated from the 
higher M reco i] region with large uncertainties. For the 
e + e~ — > D+D~ Born cross-section we set an upper limit 
of 0.04 pb at the 90% CL. 

The relative sizes of the measured cross-sections agree 
with the predictions of Ref. [l|: a(e+e~ — > D^D^~) 
and cr(e + e~ — > D+D^T) are similar, while a(e+e~ — > 
D+D~) is much smaller. e+e~ — > D*+ D*~ produc- 
tion is saturated by the D^+D*^ final state, also as 
expected. The absolute cross-sections are smaller than 
those of by a factor of 4, which is comparable to the 
theoretical uncertainty. Recent calculations based on the 
constituent quark model [3j reproduce the D* + D*~ and 
D + D*~ cross-sections very well, but predict cr(e+e~ — > 
D+D~) = 0.1 pb, somewhat larger than our limit. The 
predicted D*j+D^T fraction in e+er -> D*+D*~ produc- 
tion, 65% |3j, is smaller than we observe. 

In summary, we report the first measurement of the 
cross-sections for the e+e~ — > D^D^~ and e + e~ — > 
D+D^T processes at y/s — 10.6 GeV and set upper limits 
on the e+e~ -> D*+D*-, e+e- -> D*+D*~, e+e" -> 
D+D*£~ and e+e~ — > D+D~ cross-sections. 
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